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We use the multifractal detrended fluctuation analysis (MF-DFA) to study the electrical discharge 
current fluctuations in plasma and show that it has multifractal properties and behaves as a weak 
anti-correlated process. Comparison of the MF-DFA results for the original series with those for the 
shuffled and surrogate series shows that correlation of the fluctuations is responsible for multifractal 
nature of the electrical discharge current. 



I. INTRODUCTION 

Plasma physics is concerned with the complex inter- 
action of many charged particles with external or self- 
generated electromagnetic fields. It plays an essential 
role in many applications, ranging from advanced lighting 
devices to surface treatments for semiconductor applica- 
tions or surface layers. At the same time, the interpreta- 
tion and estimation of physical and chemical properties of 
a plasma fluid have been one of the main research areas 
in the science of magnetohydrodynamics and transport 
theory [1-15]. As in many other fields in physics, the 
complex physics requires advanced numerical tools to be 
developed and used. 

It is well-known that the discharge current fluctuations 
in the plasma often exhibit irregular and complex behav- 
ior. Indeed, the current fluctuations represent a dynam- 
ical system influenced by many factors, such as the pres- 
sure, the electrical potential between cathode and anode, 
the electrical properties of the gas, noises, and trends, 
due to the experimental setup, etc. Factors that influ- 
ence the trajectory of discharge current fluctuations have 
enormously large phase space. Thus, the use of stochas- 
tic tools for investigating their statistical properties is 
natural. Because of the complexity and stochasticity of 
the discharge fluctuations in the plasma fluid, it is gen- 
erally difficult to have access to detailed dynamics of the 
plasma ions, without paying attention to the statistical 
aspect of plasma. Therefore, there may be no remedy, 
except using stochastic analysis to investigate the evolu- 
tion and physical properties of the discharge current pro- 
duced by such fiuctuations. Also due to the limitations in 
the experimental setup for measuring the fiuctuations. as 
well as the finiteness of the available data in some cases, 
the original fluctuations may be affected by some trends 



and non-stationarities. Therefore, in order to infer valu- 
able statistical properties of the original fluctuations and 
avoid spurious detection of correlations, one must use a 
robust method which should be insensitive to any trends. 

Fluctuation of the electric and magnetic flelds of 
plasma, spectral density, logistic mapping and nonlin- 
earity of ionization wave have been investigated in [1-14]. 
Recently, Carreras et. al., have shown that the plasma 
has a multifractal nature with intermittency levels com- 
parable to the levels measured in neutral fluid turbulence 
[14]. Also, Budaev et. al, used the scaling behavior of 
structure functions and wavelet transform modulus max- 
ima (WTMM). They showed the anomalous transport of 
particles in the plasma phase attached to the turbulent 
property has multifractal nature [8,11-13]. 

Although the analysis of discharge current in the 
plasma has a long history [16-18], nevertheless, some im- 
portant issues, such as ions and electrons acceleration 
mechanism, interaction between laser and plasma, espe- 
cially from the statistical properties point of view, fractal 
features, effects of trends in small and large scales and the 
kind of correlations have remained unexplained [1-17]. 

Generally, correlated and uncorrelated time series 
could have same probability distribution function. Also 
they may have mono-fractal or multi-fractal nature. The 
mono-fractal signals can be describe by one scaling expo- 
nent. However many time series do not exhibit a simple 
monofractal scaling behavior. In some cases, there exist 
crossover (time-) scales separating regimes with different 
scaling exponents. In other cases, the scaling behavior 
is more complicated, and different scaling exponents are 
required for different parts of the series. This occurs, e. 
g., when the scaling behavior in the flrst half of the se- 
ries differs from the scaling behavior in the second half. 
In even more complicated cases, such different scaling 



behavior can be observed for many interwoven fractal 
subsets of the time scries. In this case a muhitude of 
scahng exponents is required for a fuU description of the 
scahng behavior, and a multifractal analysis must be ap- 
plied. In nature, two different types of multifractality 
in time series can be distinguished: (i) Multifractality 
due to a broad probability density function for fluctua- 
tions. In this case the multifractality cannot be removed 
by shuffling the series, (ii) Multifractality due to differ- 
ent (time-) correlations for small and large fluctuations. 
If both kinds of multifractality are present, the shuffled 
series will show weaker multifractality than the original 
series. 

Here we rely on the state-of-the-art of computational 
methods in statistical physics to characterize the complex 
behavior of electrical discharge time series. We study 
the discharge current fluctuations (see the upper panel of 
Figure 1) by the multifractal detrended fluctuation anal- 
ysis (MF-DFA) and Fourier-detrended fluctuation anal- 
ysis (F-DFA) methods. Using the method proposed in 
[19,20], we investigate the relation between the ampli- 
tude and the period of the trend and crossover in dis- 
charge current fluctuations in the plasma. 

The rest of this paper is organized as follows: In Sec- 
tion II we describe some important steps of fractal analy- 
sis to explore the stochastic time series in the presence of 
sinusoidal trends. The Hurst exponent and its relation 
to the classical multifractal, the generalized multifrac- 
tal dimension and the Holder exponents are described in 
section II. Data preparation, experimental setup, surro- 
gate and shuffled time series are also given in section II. 
We then eliminate the sinusoidal trends via the F-DFA 
technique in section HI, and investigate the multifractal 
properties of the remaining fluctuations. In section IV we 
deal with the source of multifractality in data. Section 
V is devoted to a summary of the results. 



II. ANALYSIS TECHNIQUES AND 
EXPERIMENTAL SETUP 

The simplest type of multifractal analysis is based 
upon the standard partition function, which has been de- 
veloped for multifractal characterization of normalized, 
stationary measurements [21-24]. The standard formal- 
ism does not yield the correct results for nonstationary 
time series that are affected by trends. The improved 
multifractal formalism has been developed by Muzy et 
al. [25], which is known as the wavelet-transform modu- 
lus maxima (WTMM) method [25,26]. It is based on the 
wavelet analysis and involves tracing the maxima lines 
in the continuous wavelet transform over all the scales. 
The second method which is known as MF-DFA, is based 
on the identification of the scaling behavior of the qth 
moments and is the generalization of the standard DFA 
which uses only the second moment, q = 2 [19,20,27-32]. 




FIG. 1. Upper panel: typical discharge current fluctuations 
as a function of time in our experimental setup. Lower panel: 
the sketch of the experimental setup to record the discharge 
current fluctuations in the tube, filled by helium, with (1) 
vacuum pumps; (2) copper cylinder; (3) power supply; (4) hot 
cathode; (5) anode plate; (6) glass tube; (7) single Langmuir 
probe; (8) water cooling; (9) Pirani pressure gage; (10) gas 
inlet; (11) OP-Amp; (12) low pass filter, and (13) A/D card 
and PC. 



It has successfully been applied to diverse problems, 
such as heart rate dynamics [33-35] , economical time se- 
ries [36-41], river flow [42] and sunspot fluctuations [43], 
cosmic microwave background radiations [44] , and music 
[45-47]. 

In general, experimental data are often affected by 
non-stationarities, such as trends which must be well- 
distinguished from the intrinsic fluctuations of the se- 
ries, in order to determine their correct scaling behav- 
ior. In addition, very often we do not know the rea- 
sons for the underlying trends in the collected data and, 
even worse, we do not know the scales of the underly- 
ing trends. For reliable detection of the correlations, it 
is essential to distinguish trends from the intrinsic fluc- 
tuations in data. Hurst rescaled-rangc analysis [48] and 
other non-detrending methods work well if the records 
are long and do not involve trends. But, if trends are 
present in the data, they might yield inaccurate results. 
In general, the MF-DFA is a well-established method 
for determining the scaling behavior of noisy data in 



the presence of trends without knowing their origin and 
shape [27,28,33,49,50]. Here in order to eliminate the 
effect of sinusoidal trend, we apply the Fourier DFA (F- 
DFA) [51,52]. After elimination of the trend we use the 
MF-DFA to analyze the data set. 

A. MF-DFA Method 

The MF-DFA consists of the following four steps (see 
[19,20,27,29-32,53] for more details): 

(i): Computing the profile of underlying data series, 

Xfc, as 

i 

Y{t) = Y,[xk-{x)] z = l,...,7V (1) 

fe=i 

(ii) : Dividing the profile into Ns = mt{N/s) non- 
overlapping segments of equal lengths s, and then com- 
puting the fluctuation function for each segment 

1 

F\s, m) = - y {Y[{m - l)s + i] ~ yrn{i)f (2) 

where ym{i) is a fitting polynomial in segment mth. Usu- 
ally, a linear function is selected for fitting the function. 
If there do not exist any trends in the data, a zeroth- 
order fitting function might be enough [27,28,35]. 

(iii) : Averaging the local fluctuation function over all the 
part, given by 

I m=l ) 

Generally, q can take any real value, except zero. For 
q = 2, the standard DFA procedure is retrieved. 

(iv) : The final step is determining the slope of the log- 
log plot of Fq{s) versus s directly determines the so-called 
generalized Hurst exponent h{q), as 

F,{s) ^ (4) 

For stationary time series, such as the fractional Gaus- 
sian noise (fGn), Y{i) in Eq. (1) will be a fractional 
Brownian motion (fBm), and so, < h{q = 2) < 1.0. 
The exponent h{2) is identical with the well-known Hurst 
exponent H [21,27-29]. Moreover, for a nonstationary se- 
ries, such as the fBm, Y{i) in Eq. (1) will be a sum of 
the fBm series and, thus, the corresponding scaling ex- 
ponent of Fq{s) is identified by h{q = 2) > 1.0 [27,28,54] 
(see the appendix of [42,43] for more details). In this 
case, the relation between the exponent h{2) and H is 
H = h{q = 2) — 1. The auto-correlation function is char- 
acterized by a power law, C{s) = (nkUk+s) ~ s~'^ ^ with 
7 = 2 — 2H. Its power spectrum is given by, 5'(z^) ~ , 



with frequency v and (3 = 2H — 1. In the nonstationary 
case, the correlation function is 

C{t,j) = {x,x,) ^ - 1* - jf ^ (5) 

where i,j > 1 and the power-spectrum scaling exponent 
is, /3 = 2H + 1 [27,28,42,43,54]. 

For monofractal time series, h{q) is independent of q, 
since the scaling behavior of the variances F'^{s,m) is 
identical for all the segments m, and the averaging pro- 
cedure in Eq. (3) will just yield identical scaling behavior 
for all values of q. If we consider positive values of q, the 
segments ra with large variance F'^{s,m) (i.e., large de- 
viations from the corresponding fit) will dominate the 
average Fq(s). Thus, for positive values of <j, h{q) de- 
scribes the scaling behavior of the segments with large 
fluctuations. For negative values of q, on the other hand, 
the segments m with small variance F'^{s, m) will domi- 
nate in the average Fq{s), and h{q) describes the scaling 
behavior of the segments with small fluctuations. 

The classical multifractal scaling exponents T(g), de- 
fined by the standard partition function-based formalism, 
discussed in literature [21-24,30-32,43], is related to the 
generalized hurst exponent via the MF-DFA as 

T{q) = qh{q) - 1 (6) 

Moreover, the generalized multifractal dimensions D{q) 
read as 

Diq) . ^ = (7) 
q-1 q-1 

Another way of characterizing a multifractal series is 
through its singularity spectrum /(a), which is related 
to T{q) via a Legendre transform [21,23]. Here, a is the 
singularity strength or the Holder exponent. Using Eq. 
(6), we can directly relate a and f{a) to h{q) 

a = h{q) + qh'{q) and /(a) = q[a - h{q)] + 1 (8) 

A single Holder exponent denotes monofractality, while 
in the multifractal case, different parts of the structure 
are characterized by different values of a, leading to the 
existence of the spectrum /(a). 

In some cases, there exist one or more crossover (time) 
scales, Sx, segregating regimes with different scaling ex- 
ponents, e.g., the correlation exponent for s <C Sx and 
another type of the correlation or uncorrelated behavior 
for s ^ Sx [19,20,30-32]. In the presence of different be- 
havior of the various moments in the MF-DFA method, 
distinct scaling exponents are required for different parts 
of the series [19,20]. Therefore, one needs a multitude of 
scaling exponents (multifractality) for a full description 
of the scaling behavior. A crossover usually can arise 
from a change in the correlation properties of the sig- 
nal at different time or space scales, or can often arise 
from trends in the data [19,20]. However it is well-known 



WTMM method can remove this crossover but. in many 
case, the presence of crossover as well as their values have 
physical importance. 

Let us mention two advantages of the MF-DFA method 
in order to compare with the applicability of WTMM 
method. The first advantage corresponds to the effort 
of programming and second one is related to the per- 
formance and reliability of given results. The MF-DFA 
does not require the modulus maxima procedure while 
WTMM need to do this task. The wavelet coefficients 
can become arbitrary small in WTMM method. Sub- 
sequently MF-DFA does not involve more effort in pro- 
gramming as well as more time consuming in contrast 
to WTMM, specially for long length time series such as 
our underlying data set [53]. The second advantage is 
related to the fact that the MF-DFA gives more reliable 
results than WTMM specially, for negative moments. It 
has been reported that WTMM gives an overestimated 
multifractal exponent and in some cases WTMM can also 
give different results if one applies different wavelets [55]. 
However, the most disadvantage and limitation of MF- 
DFA method will appear when it is applied to investigate 
data set with small size. In this case MF-DFA gives rich 
singularity spectrum corresponds to more multifractality 
property than what should be existed. Fortunately this 
circumstance does not occur in our situation, because the 
typical size of our plasma data is about 10®. 

In order to remove the trends correspond to the low 
frequency periodic behavior, we transform the recorded 
data to the Fourier space using the method proposed in 
[56] (see also [19,20,30-32]). Using this method we can 
track the influence of sinusoidal trends on the results and 
determine the value of so-called crossover in the fluctua- 
tion function, in terms of the scale in DFA method. We 
determine over which scale noises or trends have domi- 
nant contribution [51,52,57,58]. After removing the dom- 
inant periodic functions, such as sinusoidal trends, we 
obtain the fluctuation exponent by direct application of 
the MF-DFA. If truncation of the number of the modes 
be sufficient, the crossover due to a sinusoidal trend in 
the log-log plot of Fq{s) versus s disappear. 



B. Data description and experimental setup 

To investigate the stochastic nature of the discharge 
current fluctuations in a typical plasma, we constructed 
an experimental setup shown in the lower panel of Figure 
1. The discharge glass tube has two copper ends, 80 mm 
in diameter and 110 cm in length. One end is the anode 
electrode (a flat copper plate as a positive pole) , while the 
other end represents the cathode (tungsten fllament as a 
negative pole and electron propagator). The discharge 
tube is evacuated to a base pressure of 0.1 up to 0.8 torr 
under a voltage of 400 — 900V and flUed with Helium 
as the working gas. The pressure, voltage and current 



should be optimal for ensuring the stability of the plasma. 
The discharge current fluctuations were monitored using 
a resistor which was connected to an operational ampli- 
fier impedance converter. We fixed the pressure and in- 
vestigated how the statistical properties of plasma chang- 
ing under variation of the current. The fluctuations of 
the discharge current were digitized and cleaned with a 
filter that omitted direct current. Thereafter, the fluctu- 
ation of the discharge were recorded at a rate of 44100 
sample / sec, with a resolution of 12 bits, using a analog- 
digital card for several values of the electrical discharge 
current intensity, namely, 50,60,100,120,140,180, and 
210 mA. The typical size of the recorded data for every 
current intensity is about 10®. 



C. Surrogate and shuffled data 

The phase-randomized surrogate consists of three steps 
[59-62]: 

(i) Computing the discrete fourier transform (DFT) 
coefficients of the series 



(9) 



where i/ = k/NAt and At is the step of digitalization in 
the experimental setup. 

(ii) Multiplying the DFT coefficients of the series by 
a set of pseudo-independent, uniformly distributed 4'{h') 
quantities in the range [0, 27r) 



(10) 



(iii) The surrogate data set is given by the inverse DFT 



N-l 



fel f,i't'(k)^-i'2TTnk/N 



k=0 



(11) 



The power spectrum of the surrogate data set is the 
same as one for the original data. According to the 
Wiener-Khintchine theorem, surrogate data has the same 
autocorrelation as the initial series [59,60]. The ampli- 
tude of the surrogate data will be preserved as in the 
original data. However the probability density function 
of the data will change to the Gaussian distribution. We 
note that applying the DFT needs the data to be pe- 
riodic [63]. In addition, this procedure eliminates non- 
linearities, preserving only the linear properties of the 
underlying original data set [60] . 

To produce a shuffled data set, one should clean the 
imposed memory in series. To this end, we should ran- 
domize the order of data in underlying series while their 
values remain unchanged. 



III. FRACTAL ANALYSIS OF ELECTRICAL 
DISCHARGE FLUCTUATIONS TIME SERIES 

As mentioned in section II, spurious correlations may 
be detected if the time series be nonstationary, or is af- 
fected by trends. In such cases, direct calculation of the 
correlation exponent, the spectral density, the fractal di- 
mensions, etc., do not yield reliable results. Our data 
sets are affected by some trends, such as the alternative 
current oscillation, the noise due to the electronic instru- 
ments, and the fluctuations of striation areas near the 
anode and cathode plates. Therefore, we must use de- 
trended methods to distinguish the intrinsic fluctuations 
from the nonstationarity and trends. 

Let us determine whether the data set has a sinusoidal 
trend or not. According to the MF-DFA method, the 
generalized Hurst exponents h{q) in Eq. (4) are deter- 
mined by analyzing the log-log plots of Fq{s) versus s for 
each q. Using the rate of digitization in the experimental 
setup, 44100 sample / sec, one can simply change the unit 
of s to seconds. It must be pointed out that, to infer the 
desired exponents and to avoid errors arising for small 
values of s [19,20,30-32], we use the interval s > 0.005 
sec in our analysis. The resolution of the recorded data 
in our setup is 1/44100 ~ 0.00023 sec. We use this inter- 
val throughout the paper, unless specified otherwise. Our 
investigation indicates that there is at least one crossover 
time scale in the log- log plots of Fq{s) versus s for ev- 
ery q. To determine its value, we use the following two 
criteria and combine their results: 

(i) Based on the recent results by [19,20] and [42,43], 
every sinusoidal trend in the data causes some crossovers 
in the scaling function, Fq{s), derived by the MF-DFA. 
The number of such crossovers depends on the size of 
data and the wavelength of the sinusoidal trends [19,20]. 
It is well known that the crossovers divide the fluctua- 
tions function into some regions that correspond to var- 
ious scaling behaviors of Fq{s) versus s, that are related 
to the competition between noise and trends [19,20]. To 
prove this statement and show how one can determine the 
value of crossovers, we generated numerically a time sc- 
ries which is a superposition of a correlated noise, namely, 
with the Hurst exponent H = 0.8, and a sinusoidal trend 
with its period equal to T = 20 sec. As shown in Fig- 
ure 2, one sinusoidal trend is embedded in the data, and 
there are two crossovers at S3 ^ 3 sec and S4 ~ 20 sec in 
the fluctuation function. The larger crossover is directly 
associated with the period of the sinusoidal trend [19,20]. 

The crossovers are also confirmed by resorting to the 
power spectrum of the data shown in the middle graph 
of the lower panel of Figure 2. One can observe that for 
1^ > ~ I/S3 ~ 0.3 scc~^, the scaling behavior of the 
power spectrum is the same as those for the correlated 
noise, indicating that in these scales the noise effect is 
dominant. Consequently, the scaling behavior of Fq{s) 



for s < S3 is very close to those of the correlated noise 
alone (see the upper panel of Figure 2). 

Thereafter, we embedded ten sinusoidal trends with 
various frequencies in the original signal, and performed 
the same computations as shown in Figure 2. We found 
that, when one increases the number of sinusoidal trends 
in the original noise, then expects that the value of the 
crossover at large scales, namely, s^lvi) ( which is related 
to the dominant embedded sinusoidal trends), will ex- 
tend, as shown in Figure 2. In other words, we have an in- 
terval, i.e., s e [s2, S4], or in frequency space, v G [vi, 1/3], 
within which the scaling behavior of the fluctuation func- 
tion changes smoothly. Therefore, we cannot determine 
an exact value for the crossover at large scales. More- 
over, as shown in the upper panel of Figure 2, at time 
scales smaller than S2 or for frequencies larger than 1/3, 
the fluctuation function retrieves its noisy behavior. This 
observation depends on the amplitudes and frequencies 
of the embedded periodic trends. Therefore, the exact 
value of the crossover at small scale is not obvious in the 
spectral density and, hence, we use the following criteria. 

ii) As discussed above the existence of many sinusoidal 
trends in our data set, we expect that a plot of Fq{s) 
versus s possess at least one crossover. This crossover 
divides Fq{s) into two regions, as shown in the upper 
panel of Figure 3 (for instance, we take q = 2 and use 
the data set with 7 = 50 mA) . To determine the value of 
the crossover, we introduce a A(s) function as: 

A(s) = ^/[F\^i^^F[~^ (12) 

for each g, where F{s) and i^Lincar(s) are the fluctuation 
functions for the original data and the filtered data pro- 
duced by the F-DFA method (see below), respectively. 
In Figure 4, we plot A(s) as a function of s for the 
plasma fluctuations. Hereafter, for convenience we omit 
the subscript q and take q = 2 unless expressed other- 
wise. The crossover occurs at Sx ~ 0.02 sec, correspond- 
ing to Vx ~ 50 sec~^. Clearly, the fluctuation function 
for s < Sx has the same scaling behavior as the noise 
without trends (see Figure 3). 

As expressed in the last section, to cancel the sinu- 
soidal trend in the MF-DFA, we have applied the F-DFA 
method. Indeed, we truncate the lowest frequencies up 
to the value that the regression of the linear fitting of 
the corresponding log-log plot of the fluctuation function 
for truncated data reaches 0.95. It is well known that if 
we truncate the frequencies more than the necessary val- 
ues for eliminating the crossover, some statistical prop- 
erties of the underlying noise might be lost [19,20,52]. 
To eliminate the crossover scales, we need to remove ap- 
proximately the first 400 terms of the Fourier expansion. 
Then, by the inverse Fourier transformation, the fluctu- 
ations without the sinusoidal trend are retrieved. The 
result is shown in Figure 5. The generalized Hurst expo- 
nent, the classical multifractal scaling exponents, and the 
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FIG. 2. Upper panel shows the fluctuations function for 
the correlated noise and its superposition with one and ten 
sinusoidal trends. Lower panel corresponds to their power 
spectrum. The curves, from top to bottom, correspond to the 
correlated noise, one and ten embedded sinusoidal trends in 
the noise, respectively. Oblique solid lines correspond to the 
scaling behavior of power spectrum of the clean correlated 
noise, namely, [3 = 0.6 in S{u) 
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singularity spectrum for the data after the elimination of 
sinusoidal trends are illustrated in Figure 6. In addi- 
tion, we estimated the errors at la confidence level for 
all the derived values reported in all the tables, using the 
common statistical method [64] . The cleaned time series 
is a multifractal process, as indicated by the strong q- 
dependence of the generalized Hurst exponents [53] . The 
g-dependence of the classical multifractal scaling expo- 
nent T{q) has different behaviors for q < and q > 0. 
For both positive and negative values of q, the value of 
slopes of T{q) are indicated in Figure 6. As mentioned 
before for q < 0, small fluctuations will be dominant in 
the fluctuation function, whereas for g > the large fluc- 
tuations represent the dominant effect in Fq{s). From the 
statistical point of view, usually, for a multifractal anti- 
correlated series, namely, one with H < 0.5, the value 
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FIG. 3. Upper panel shows the crossover behavior of the 
log-log plot of F{s) versus s for the original time series for 
q = 2.0. Lower panel corresponds to the power spectrum of 
the original and cleaned data set. Oblique solid lines corre- 
spond to the scaling behavior of the power spectrum of the 
cleaned series, namely, (3 = —0.16 in S{u) ~ . 



of h{q) for (7 < is smaller than the generalized Hurst 
exponent for positive moments. This is due to the fact 
that, the number of the large fluctuations are statisti- 
cally much more than the small fluctuations in the time 
series. In other words, the series is intermittent. The 
same circumstances arise for the plasma fluctuations. In 
the presence of free charge, every large deviations in the 
electrostatic equilibrium would be shielded by a cloud 
of oppositely charged particles [49,65,66]. Therefore, we 
expect the plasma fluctuations to be an anti-correlated 
series. The value of the Hurst exponent confirms that 
the data set is a stationary process. According to the 
MF-DFA results, all of the discharge current series be- 
have as weak anti-correlated process. The value of the 
Hurst exponent and the classical multifractal exponent 
for g = 2, in the region that the sinusoidal trend is not 
pronounced, calculated via the MF-DFA method are re- 
ported in Table I. The correlation and power spectrum 
exponents are given in Table II. 




FIG. 4. The function A(s) for the plasma fluctuations ver- 
sus s. 

In spite of the power of the MF-DFA method, some 
cases encounter problems and the method yields inac- 
curate results. The DFA method can only determine 
a positive Hurst exponent, but yields an inaccurate re- 
sult for the strongly anti-corrclatcd recorded data, when 
if ~ 0. To avoid this situation, one should use the inte- 
grated data. In that case the scries is the so-called double 
profiled data set. The corresponding Hurst exponent is 
H = H — 1, where H is derived from the DFA method 
for the double profiled series [43,53]. 

Let us now discuss the finite size effect of the data 
set on the results. According to the recent analysis in 
Ref. [30-32], a deviation from the DFA results occurs 
in short records. The modified version of the MF-DFA 
should be used for such cases [30-32]. Usually, in the 
MF-DFA method, the deviation from a straight line in 
the log-log plot of Eq. (4) occurs for small scales s. The 
deviation limits the capability of DFA for determining the 
correct correlation behavior at very short scales, and in 
the regime of small s. The modified MF-DFA is defined 
as follows [30-32] 



71 mod 
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[F«huf(s)]2)l/2s'l/2 



(13) 



where ([-F|^"''(s)]^)^/^ denotes the usual MF-DFA fiuctu- 
ation function [defined in Eq. (3)] , averaged over several 
configurations of the shuffled data taken from the original 
time series, and s' « A^/40. The values of the Hurst expo- 
nent obtained by the modified MF-DFA methods for the 
time series are reported in Table HI. The maximum rela- 
tive deviation of the Hurst exponent, which is computed 
by the modified MF-DFA, relative to the MF-DFA for 
the original data, is approximately 2.7%. Moreover, Fig- 
ure 7 shows a comparison between the generalized Hurst 
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FIG. 5. The MF-DFA functions F,(s) in terms of the time 
scale s, in a log-log plot. The original time series, m = 0, 
the truncation of the first 200, m = 200, and 400 terms, 
m — 400. This plot is for a typical value of the electrical 
discharge current intensity, / = 50 mA. 



exponent, derived by the common MD-DFAl and modi- 
fied MF-DFA. It indicates that the modified MF-DFA is 
consistent with the MF-DFA at Icr confidence level for 
various moments. 

By inspecting the log-log plot of the fiuctuation func- 
tion versus s (i.e., F2{s) = C-hs^), we find the depen- 
dence of its amplitude on the electrical discharge current. 
Indeed, Rcf. [42] showed that this amplitude is given by 

^2 4^2 



— 

^ {2H+1) 2H + 2 
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H + 1 2H + 1 
3^2 



1 - 



'H + l){2H+l] 



(14) 



where = (.T(i)2^ is the variance of the data. Figure 
8 shows the value of the amplitude as a function of the 
current. 



IV. SEARCH FOR ORIGIN OF 
MULTIFRACTALITY IN THE DATA SET 

Due to the strong q dependence of the generalized 
Hurst exponent and the distinct slopes T{q) for the vari- 
ous moments, the remaining data set, after the elimina- 
tion of the sinusoidal trends, has multifractal properties. 
In this section we search for the source of the multifrac- 
tality. In general, two different types of multifractality 
in time series may exist: 

(1) Multifractality due to the fatness of the probability 
density function (PDF) of the time series, in compari- 
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FIG. 6. Upper panel shows the generahzed Hurst exponent 
versus q for the original, the surrogate, and the shuffled se- 
ries without the sinusoidal trend. Lower panel left and right 
correspond to the classical multifractal scaling exponent and 
the singularity spectrum f{a) for the data set at 7 = 50 niA, 
respectively. 



son to a Gaussian PDF. In this case, the multifractality 
cannot be removed by shuffling the series, because the 
correlations in the data set are affected by the shuffling, 
vifhile the PDF of the series is invariant. 

(2) MultifractaUty due to different types of correlations 
in the small and large scale fluctuations. In this case, 
the data may have a PDF with finite moments, e.g., a 
Gaussian distribution. Thus, the corresponding shuffled 
time series will exhibit monofractal scaling, since all the 
correlations arc destroyed by the shuffling procedure. If 
both kinds of multifractality are present, the shuffled sc- 
ries will exhibit weaker multifractality than the original 
series. 

The easiest way to distinguish the type of multifractal- 
ity is by analyzing the corresponding shuffled and surro- 
gate time series. The shuffling of the time series destroys 
the correlation. Therefore, if the multifractality only be- 
longs to the correlation, we should find /ishuf(<z) = 0.5. 
The multifractality due to the fatness of the PDF series 
is not affected by the shuffling procedure. On the other 
hand, to determine whether the multifractality is due to 
the broadness or fatness of the PDF, the surrogate data 
are used [67]. If multifractality in the time series is due 
to a broad PDF, h{q) for surrogate data will be indepen- 
dent of q. If both kinds of multifractality are present in 
fluctuations, the shuffled and surrogate series will exhibit 
weaker multifractality than the original one. The utility 
of the two tests was first demonstrated by Ivanov et al. 
[68-70]. 

To check the nature of the multifractality, the we com- 
pare the fluctuation function Fq{s) for the original series 
(after removal of the sinusoidal trends) with the results 




FIG. 7. Comparison between the modified and commonly 
used MF-DFA for a typical value of the electrical-discharge 
current intensity, / = 50 mA. 

for the corresponding shuffled, i^|^"^(s) and that of the 
surrogate series, F™''(s). The differences between the 
two fluctuation functions with the original one, directly 
indicate the presence of the correlations or broadness of 
the PDF in the original series. The differences can be ob- 
served in plots of F,(s)/F|'^"f(s) and F,(s)/F|"(s) ver- 
sus s [53]. Since the anomalous scaling due to a broad 
PDF affects Fq{s) and F^^^^{s) in the same way, only 
multifractality due to the correlations will be observed 
in Fq{s) / Fq^^^ (s) . The scaling behavior of the two ratios 
arc given by 
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If only the fatness of the PDF is responsible for 
the multifractality, one should obtain, h(q) = /ishuf('7), 
and, hcor{<l) = 0. On the other hand, deviations from 
hcoril) ~ indicates the presence of correlations, while 
the q dependence of /icor('z) indicates that multifractality 
is due to the correlation. If both the distribution and 
correlation multifractality are present, both /ishuf (<?) and 
/isur(<?) will depend on q. 

The q dependence of the exponent h{q) for the original, 
surrogate, and shuffled time series are shown in Figures 6. 
The q dependence of hpuF shows that the multifractal- 
ity nature of the time series is almost due to the correla- 
tion. However, the value of /ipdf('z) is deviated from zero 
which confirms the multifractality due to the broadness 
of the PDF is much weaker than the multifractality due 
to the correlation. The deviation of hsuril) and /ishuf('?) 
from h(q) may be determined by using the test: 
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FIG. 8. The value of Cn as a function of the discharge 
current intensity. 



where the symbol o can be replaced by " sur" and " shuf ' , 
in order to determine the confidence level of hgur and 
hshui, the generalized Hurst exponents of the original se- 
ries, respectively. The value of the reduced chi-square, 
xi = X^/A/" where Af is the number of degree of free- 
dom, for the shuffled and surrogate time scries arc shown 
in the upper panel of Figure 9. On the other hand, the 
lower panel of Figure 9 illustrates the width of the sin- 
gularity spectrum, /(a), i.e., Aa = a{q,nin) - a{qmax), 
for the original, surrogate, and shuffled data set. These 
values also confirm that the multifractality due to the 
correlations is dominant [55]. 

The values of the generalized Hurst exponent h{q = 
2.0), and the multifractal scaling T{q = 2) for the origi- 
nal, shuffled, and surrogate of the discharge ffuctuation 
obtained with the MF-DFA method are reported in Table 
L 



V. CONCLUSION 

Discharge current ffuctuations in plasma are affected 
by many factors. From a statistical point of view, and 
in order to understand the complexity of the ffuctua- 
tions, we applied a robust method, such as the dctrend- 
ing fractal analysis to infer the complexity and multi- 
fractal features of the underlying plasma fluctuations. In 
the presence of nonstationarity, non-detrending methods 
will be encountered with some challenges, such that they 
yield wrong or at least inaccurate results. Multifractal 
detrended fluctuation analysis is well established for in- 
vestigating noisy time series, and can be used to gain 
deeper insight into the processes that occur in nonsta- 
tionary dynamical systems, such as electrical discharge 
current. 



TABLE I. Values of H = h{q = 2) and the classical mul- 
tifractal scaling exponents for q — 2.0 for the original, sur- 
rogate and shuffled data set and different electrical currents, 
obtained by the MF-DFA. 



Sample 




H 


r 




Original 


0.42 ±0.01 


-0.16 ±0.02 


50mA 


Surrogate 


0.40 ± 0.01 


-0.21 ± 0.02 




bnutned 


n dpi 1 n n 1 

0.50 ± 0.01 


0.004 ± 0.020 




Original 


0.45 ±0.01 


-0.10 ±0.02 


60mA 


Surrogate 


0.42 ± 0.01 


-0.16 ± 0.02 




bhutned 


n dpi 1 n n 1 

0.50 ± 0.01 


0.002 ± 0.020 




Original 


0.37 ±0.01 


-0.25 ± 0.02 


100mA 


Surrogate 


36 ± 01 


—0.28 ± 0.02 




Shuffled 


0.49 ±0.01 


0.00 ±0.02 




Original 


0.38 ±0.01 


-0.23 ±0.02 


120mA 


Surrogate 


0.36 ±0.01 


-0.28 ±0.02 




Shuffled 


0.50 ±0.01 


0.00 ±0.02 




Original 


0.41 ±0.01 


-0.17 ±0.02 


140mA 


Surrogate 


0.40 ±0.01 


-0.21 ±0.02 




Shuffled 


0.50 ±0.01 


0.00 ±0.02 




Original 


0.45 ±0.01 


-0.09 ±0.02 


180mA 


Surrogate 


0.40 ±0.01 


-0.20 ±0.02 




Shuffled 


0.50 ±0.01 


-0.003 ±0.020 




Original 


0.48 ±0.01 


-0.04 ±0.02 


210mA 


Surrogate 


0.44 ±0.01 


-0.13 ±0.02 




Shuffled 


0.51 ±0.01 


0.01 ±0.02 



TABLE II. Values of the correlation and power spectrum 
exponents for the original data set with different electrical 
currents, obtained by the MF-DFA. 



Sample 


7 


P 


50mA 


1.16 ±0.02 


-0.16 ±0.02 


60mA 


1.10 ±0.02 


-0.10 ±0.02 


100mA 


1.26 ±0.02 


-0.26 ±0.02 


120mA 


1.24 ±0.02 


-0.24 ±0.02 


140mA 


1.18 ±0.02 


-0.18 ±0.02 


180mA 


0.10 ±0.02 


-0.10 ±0.02 


210mA 


1.04 ±0.02 


-0.04 ±0.02 



TABLE III. Values of Hurst exponent using the regular 
MF-DFA and modified MF-DFA for data sets. 



Sample 


H 




50mA 


0.42 ±0.01 


0.42 ±0.01 


60mA 


0.45 ±0.01 


0.45 ±0.01 


100mA 


0.37 ±0.01 


0.38 ±0.01 


120mA 


0.38 ±0.01 


0.37 ±0.01 


140mA 


0.41 ±0.01 


0.42 ±0.01 


180mA 


0.45 ±0.01 


0.46 ±0.01 


210mA 


0.48 ±0.01 


0.48 ±0.01 
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FIG. 9. Upper panel shows to the values of xt for the data 
sets. Lower panel shows the width of the singularity strength, 
Aa, for the original, surrogate and shuffled data sets in vari- 
ous electrical discharge currents, obtained by the MF-DFA. 



We showed that the MF-DFA resuhs for the time se- 
ries for various electrical currents have one crossover time 
scale, Sx- This crossover time scale is about Sx ~ 0.02 
second, and is related to the sinusoidal trends. The 
crossover time scale which discriminates the noise and 
trends intervals can be potentially related to the coher- 
ent time scale in turbulent plasma. Plasma fluctuations 
are not affected by external factors within time scale. 
To minimize the effect of trends and produce clean data 
set for further investigation, we applied the Fourier de- 
trended fluctuation analysis to the data sets. Indeed, af- 
ter applying the F-DFA, the data set without sinusoidal 
trends is recovered, and the spurious behavior in the MD- 
DFA results disappear. Applying the MF-DFA method 
on the cleaned data set demonstrated that the discharge 
current fluctuations are stationary time series. 



According to the value of the Hurst exponent, com- 
puted by the MF-DFA method, all the discharge cur- 
rent time series behave as weak anti-persistent processes. 
These findings can be interpreted as follows: in the pres- 
ence of free charges, every large deviation from the elec- 
trostatic equilibrium is shielded by a cloud of oppositely- 
charged particles [49,65,66]. This also may be related to 
the fast dissipation of turbulent kinetic energy in plasma. 
Our results also confirmed that the multifractality nature 
is a global property of various plasma data based on dif- 
ferent experimental setups [11]. We found that Hurst 
exponent and multifractality nature based on singularity 
spectrum didn't depend on the discharge current inten- 
sity (see Table HI and Figure 9). This result indicates 
that increasing the amount of charged particles at least 
in our experimental setup almost don't alter the statisti- 
cal properties of the plasma fluid. But it is interesting to 
extend these analysis to a broader set of plasma data in 
various working pressure and check their statistical prop- 
erties. The (/-dependence of h{q) and T{q) indicated that 
the data sets have multifractal properties. 

The value of h{q) for g > is larger than the same 
quantity for g < 0, indicates that the number of large 
fluctuations are statistically larger than the small fluctua- 
tions in the time series. Our results show that, amplitude 
of the fluctuation function Cn is a monotonous function 
in terms of discharge current intensity. This demonstrate 
that by increasing the current intensity, plasma instabil- 
ity will be occurred [13] and consequently one has large 
variance for data set (see Figure 8). 

In order to recognize the nature of the multifractality, 
we compared the generalized Hurst exponent of the orig- 
inal time series with those of the shuffled and surrogate 
ones. The comparison indicated that the multifractal- 
ity due to the correlations makes more significant con- 
tribution than the broadness of the probability density 
function of the current fluctuations. 
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